Dependence of asymmetries for charge distribution with respect 
to the reaction plane on initial energy in heavy ion collisions 

V.A. OkorokcnQ 

National Research Nuclear University "MEPhI", Kashirskoe Shosse 31, 115409 Moscow, Russia 

(Dated: January 8, 2013) 

In the paper two combinations of correlators are denned in order to investigate the evolution of 
possible V/CV invariance violation in strong interactions with initial energy for heavy ion collisions. 
These combinations correspond to absolute and relative asymmetry of distribution of electrically 
charge particles with respect to the reaction plane in heavy ion collisions. Energy dependence of 
the parameters was derived from data of STAR and ALICE experiments. Significant decreasing 
both absolute and relative asymmetry is observed at energies ^/ijvjv < 20 GeV. This feature agrees 
qualitatively with other results of stage-I beam energy scan program in STAR experiment. General 
ff) ■ behavior of dependence of absolute asymmetry on initial energy agree reasonably with behavior of 

t-H ' similar dependence of Chern - Simons diffusion rate calculated at different values of external Abelian 

magnetic field. The observed behavior of understudying parameters vs energy can be considered 
04 ■ as indication on possible transition to predominance of hadronic states over quark-gluon degrees of 

(-H ' freedom in the mixed phase created in heavy ion collisions at intermediate energies. 
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Vacuum of quantum chromodynamics (QCD) as non- Abelian gauge theory contains topologically non-trivial con- 
figurations of gauge fields which deeply relate with V /CV invariance of strong interactions. The non-trivial topology 
of QCD vacuum opens the possibility for existence of metastable domains which possess of various properties with 
i— respect to the discrete V/CV symmetries. According to predictions of the theory at finite temperature [l| decays of 
such domains or classical transitions (sphalerons) between them in deconfinement phase of color charges with restored 
chiral symmetry can results in local topologically induced violation of V /CV invariance in strong interactions - VTTP 
effect. The one of possible mechanisms for manifestation of such domains and local TXV violation experimentally 
is azimuthal asymmetry of distribution of electric charges in final state driven by a external Abelian magnetic field 
- the so called chiral magnetic effect - CM.E Q. Experimentally charge-separation effect in relativistic heavy ion 
collisions was observed by STAR collaboration at RHIC for the first time [|| and was confirmed in further by ALICE 
* ' collaboration at LHC [ij . As indicated above deconfinemet state of color charges and restored chiral symmetry are es- 
sential for the experimental manifestation of possible local strong V/CV violation, i.e. for CMS. The former is needed 
to separate (anti-)quarks with opposite electric charges by a distance larger than nucleon size. The restored chiral 
symmetry is required because charge separation is possible at conserved chirality only. Therefore the investigation 
of energy dependence of CMS via definition of convenient parameters and analysis of its corresponding dependence 
is important challenge for study of strongly interacting matter under extreme conditions and search for transition 
domain from prevalence of hadronic states to domination of deconfinement phase of quarks and gluons. 



II. METHOD AND VARIABLES 



With taking into account possible local strong V /CV violation the invariant distribution of final state particles with 
certain sign of electrical charge a (a = +, — ) can be written as following 



E 



d 3 N a 1 d 2 N a 
dp 2-7T p±dp±dy 



1 + 2 {v nta cos (nA(f>) + a n , a sin (nA</>)} 



A</> = </> - 1< RP . 



Here cj> is an azimuthal angle of particle under study, ^rp - azimuthal angle of reaction plane, v n , a - collective flow of 
7i-th order, the parameters a n a describe the effect of V/CV violation. The expectation value of a n , a will vanish if the 
one-particle distributions are averaged over event sample because the distributions are averaged over many domains 
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and direction of charge separation changes event by event due to random sign of the Chern - Simons number (Nqs) 
of the local domain 0,0]. Therefore it would be emphasized that the CM.E is the collective effect and its investigation 
is possible via correlation analysis only. According to the theory, the correlator contained the contribution of possible 
V/CV violation effect only is given by (K^ Q/3 ) = (a nta a n ,p), where a, (3 - the electric charge signs of secondary 
particles. Only the first harmonic coefficient was analyzed so far because as expected the a\ t0 accounts for most of 
the effect under study Q. Thus the following notation (K.J a g) = (K^) is used below. Experimental correlator has 
been proposed in Q and defined as following (K?g) = (cos (4> a + (f>p — 2^^)). This observable is sensitive to the 
effect of possible local strong V/CV violation and measures the charge separation with respect to the reaction plane. 
Both theoretical and experimental correlators are averaged over pairs of particles under study in event and over all 
events in sample. It should be mentioned that correlators defined above are V-even quantities therefore the (K?g) 
may contains contributions from background effect unrelated to possible local strong V/CV violation. Theoretical 
and experimental correlators are related by following equation 0, a, [1] 

(K^) = [<«i,„«i, /J )+A$]-<K^). (1) 

Here v\ iCC is the directed flow parameter, A<f> = $; n — <&out _ the difference between background contribution of the 
in-plane correlations (3>i n ) and background contribution of the out-of-plane correlations ($ ut)- As a general rule 
the experimental contribution given by the first term on r.h.s. of (|TJ) be neglected because directed flow averages 
to zero at rapidities symmetric with respect to midrapidity region which is considered in most relativistic heavy ion 
experiments [3|, Q. I n the framework of CM.E model Q the theoretical correlator is given by 

< K -> = rot = mk> (f / /™->i 2 ; m 

f s ± n 

Here e±± = -0.5e ±T = 0.5, £±±{x±) = £) i=+i _ £± T (%x) = H i=+ _ Ci(^x), «s ~ renormalized strong 

coupling constant, N a - multiplicity of particles with electric charge sign a in event, R - radius of colliding nuclei, 
qf - electric charge (in units of e) of quark with flavor /, A Q is the difference between number in fixed charge a on 
each side of the reaction plane, numerical coefficient re ~ 1, xj_ - position of nucleon from one of the colliding nuclei 
in the plane perpendicular to the beam axis, S± — area of overlap region of colliding nuclei in the transverse plane, r 
- proper time, B - strength of external Abelian magnetic field, functions £±(x±) take into account the suppression 
effect for transitions with change of Chern - Simons number in various parts of the volume of created hot matter 
@, El • Analytic approaches for ([2]) corresponded to various scenarios of initial time and comparison with STAR results 
were discussed in @ . The most significant discrepancy between experimental and phenomenological results observed 
for semi-central collisions @ requests the additional quantitative study and clear explanation. The STAR results for 
azimuthal correlations of high p t particles with respect to the reaction plane 0] show that difference between such 
correlations is largest for semi-central collisions. Possibly, similar conclusion is valid for particles with moderate px 
too. Thus at qualitative level one can suggests that first term part A$ in r.h.s. of ([TJ can influences on (K^g) and 
results in the disagreement between experimental and phenomenological data points. 

In the framework of CA4£ model for ideal case of chiral limit and an extremely large magnetic field the average 
square of charge difference (in units of e) between opposite sides a reaction plane can be estimated as following 



(Q 2 )-4(A 



(3) 



where A t is the difference between the number of raising and lowering of Chern - Simons number transition 0. 
Let ATVcg be a total change of Chern - Simons number, which is dominated by the rate of corresponding change. 
In this paper it has been suggested that at finite temperature (T) the transition rate between vacuum states with 
different Nqq, is dominated by transition rate due to sphalerons only, without consideration of exponentially suppressed 
transitions due to instantons. But it would be mentioned that some additional study and justification may be required 
for this approach, especially in temperature range T c < T < 3T C [8j , where T c — temperature of phase transition to the 
deconfinement state of color charges. It seems temperatures reached experimentally in relativistic heavy ion collisions 
are in the range indicated above for present and even future colliders. Thus it is suggested that the rate of change of 
Nqs is the classical transition rate, i.e. the Chern - Simons diffusion rate (Fes)- O ne can suggest A t — AiV^g because 
the most probable transitions are those which connect two neighboring vacua and transitions with simultaneously 
changing of absolute value of Chern - Simons number larger than on unit are suppressed ■ By definition the T^g 
quality is the probability of Chern - Simons number changing process to occur per unit 3-dimensional volume and 
per unit time: T^g = liniy-j-cx) V _1 ([AA^Qg] 2 ), where V is the 4-dimensional volume. The parameter T^g has been 
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intensively studied in the framework of electroweak theory. The Chern - Simons diffusion rate was computed for weak 
coupling domain. The holographic AdS/CFT correspondence [|| makes it possible to derive this parameter in Af = 4 
supersymmetric Yang - Mills (SYM) plasma for limit of large number of colors (N c — > oo) in the strong coupling 
regime [10| . Then with taking into account the sphaleron approach the Chern - Simons diffusion rate in absent of 
external Abelian magnetic field can be given by 

r C s - fix)? , fix) - | x 2 /(256 ^ 3); at x L 

Here x = 9 2 N C is the 't Hooft effective coupling, g - gauge constant, the x ^ 1 corresponds to the weak coupling 
domain and x 3> 1 - to the strong coupling regime fiol ] . The latter corresponds to the relativistic heavy ion collisions 
because the created matter seems strongly correlated system. Numerical (large) coefficient k is of nonperturbative 
nature and can be found by a lattice simulation. Estimations demonstrate that the k ~ 10 for SU(2) [ll[ and k ~ 10 2 
for SU(3) [3, i.e. k ~ lO^ 1 for SU(AT) with N = 2,3. Thus r cs oc T 4 . In accordance with definition of T cs and 
([3]) one can obtain (Q 2 ) oc T^g. One the other hand in the most ideal case under study, i.e. for contribution of the 
lowest Landau level only, one can derive 



|(A±A±) T (A±A T )| = ^([AiV cs ] 2 ), <p = (2V)" 



d 4 x[^(x ± )±£ + (x ± )} 2 . 



The choice of difference of correlators in l.h.s. seems reasonable in order to corresponding quality be positive, then 

(K±±> " <K±=f) « ([AiV cs ] 2 ) cx r cs . (4) 
Taking into account the ([1]) one can define the following quantity 

A a = -[(Kl ± )-{Kl T )}, (5) 

which is absolute asymmetry for azimuthal distribution of electric charges in final state and can be suggested for 
experimental estimation of ([AiVpg] 2 ). The relative asymmetry parameter for azimuthal distribution of electric 
charges with respect to the reaction plane in event is defined as following 

_ |(Kg ± )|-|(K^)| 

|(Kl ± )| + |(K^)|- W 



III. RESULTS 



Fig. Q] shows the energy dependence of A a for semi-central events in various bins of centrality for heavy ion beams. 
Collisions of moderate Cu nuclei are characterized by significantly larger values of A a than that for heavy ions (Au, U) 
at corresponding energies. It would be noted that for Au+Au at ^snn — 62.4 and 200 GeV absolute asymmetry 
([5]) was calculated based on the preliminary STAR correlators for increased (full) statistics at former energy and for 
Run 7 data sample at second initial energy [l4| . These preliminary STAR results allow to decrease the divergence of 
experimental points and get slightly smooth dependence A a on ^/sjvjv- The equation for Chern - Simons diffusion rate 
for case of finite B was derived in [15[ . Background magnetic field typically created in relativistic heavy ion collisions 
is characterized by the strength B ~ T 2 [IfS [l^]. The regime of finite external magnetic field is investigated in 
the paper and it is observed that the curves for B = and for case B ~ T 2 coincide completely for all centrality 
bins under consideration, i.e. for full centrality domain 20-60% shown at Fig. [TJ But relations (|4]) were obtained for 
extremely strong magnetic field. Thus energy dependence for the Chern - Simons diffusion rate are computed for 
external Abelian magnetic field with various strengths. Analytical function suggested in [TtJ for description of energy 
dependence of chemical freeze-out temperature agree with available experimental data quite reasonable for energies 
up to y /s n n = 200 Gev [ijj and for all centralities under study. Therefore the analytic dependence T(^/snn) 
from [ITj is used for estimations of T^g in energy domain under consideration. The smooth curves at Fig. Q]are the 
energy dependence of Chern - Simons diffusion rate in the strong coupling regime at B — (solid), B — 5T 2 (dashed) 
and B = 10T 2 (dotted). The norm for solid curves is the STAR point for Au+Au at ^/s]^ = 200 GeV. The shaded 
bands for Chern - Simons diffusion rate at B = are defined by uncertainties of T value at fixed initial energy due 
to errors of parameters in analytic function described of T(^snn) experimental dependence [l7| . 

As seen, phenomenological curves T^s (-Jsnn) are m area caused by spread of T values at fixed ^snn f° r following 
wide range of changing of strength of external Abelian magnetic field B < 5T 2 at any initial energies under study, 
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moreover at ^/sWn 12 GeV - for all range of changing of B under considered on Fig. [TJ Functional behavior of 
^Cs(\/ s nn) does not depend on B in high energy domain at changing of B on order of magnitude at least. On the other 
hand the Chern - Simons diffusion rate decreasing faster at ^snn < 30 GeV for larger strength of external Abelian 
magnetic field. Possibly, the equation for Chern - Simons diffusion rate without external B is valid for any initial 
energies up to ^/snn = 2.76 TeV and for r > 0.1 fm/c due to uncertainty area for Fqs(^/snn) at B = discussed 
above and B values reached in heavy ion collisions [6j. Therefore further decreasing of parameter uncertainties for 
T(^/snn) behavior is essential for qualitative study of T^g dependence on strength of external Abelian magnetic field 
in nucleus-nucleus collisions. 

One sees the some decreasing of A a for LHC energy as compared to initial energy domain ^snn ~ 100 GeV, the 
absolute value of decreasing growths for more peripheral events. Absolute asymmetry of azimuthal distribution of 
electric charges in final state A a goes down significantly with ^snn decreasing for intermediate initial energy domain 
7.7-20 GeV in semi-central heavy ion collisions (Fig. [T]). For each centrality bin under study the preliminary expe- 
rimental results for A a (\/sjv]v) correspond to the energy dependence of r^g at B = on qualitative level with some 
enhancement of experimental points over phenomenological curves in range ^snn ~ 20 — 40 GeV for events with 
20-30% (Fig. [TJi) and 30-40% (Fig. [Tb) centrality. Perhaps, for better description of A a (^/sNN) at y/s^N ~ 20 - 40 
GeV the influence of external B on T^g should be taken into account. This suggestion agrees with dependence of B 
on initial energy Q which demonstrates that B in intermediate energy range ^/snn ~ 20 — 40 GeV reaches very large 
values with respect to the strength of B at ^snn ~ 100 GeV. The qualitative agreement is observed between values 
of parameter ([5]) calculated via preliminary STAR results for Au+Au collisions and phenomenological curves for any 
B under consideration. This feature can be considered as some evidence of validity of equations §5§ and dU obtained 
at ultimate strong external Abelian magnetic field for heavy ion collisions. Decreasing of ^cs(\/ s nn) observed at 
^/sjvtv < 20 GeV should leads to attenuation of CMS and its manifestation on experiment at intermediate energies. 
Taking into account conditions essential for CMS one can supposes the following hypothesis. The changing of behavior 
of dependence ^4 a (v / ^v r /v) observed at transition from high energy domain down to intermediate energy range may 
be driven by predominance of hadronic colorless states over quark-gluon deconfinement phase at y / sjvw < 19.6 GeV 
and, as consequence, by decreasing of CMS. Thus behavior of experimental parameter ([5]) vs collision energy agrees 
with qualitative expectation for transition to the predominance of hadronic phase in domain ^snn < 19.6 GeV and 
with decreasing of Chern - Simons diffusion rate at intermediate initial energies. 

Fig. [5] demonstrates the energy dependence of A r for semi-central (20-60%) events for various ion beam collisions. 
As well for Fig. [T] the dependence under consideration is some improved by calculations of parameter ^ based on 
increased STAR statistics at ^/snn = 62.4 and 200 GeV. One sees the discrepancy of A r values for moderate and 
heavy nuclei at same energies increases for more peripheral events. Behavior of A r (y^ny) f° r heavy ions changes 
dramatically at ^snn < 19.6 GeV, especially for last three bins of centrality (Fig. [5b-d). The sharp decreasing 
of relative asymmetry ^ for energies indicated above can be considered, possibly, as evidence for (significant) am- 
plification of hadronic phase influence on experimental observables. The observed feature of v4 r ( v /s]vjv) and above 
suggestio n ag ree with other STAR results obtained in the framework of stage-I beam energy scan (BES-I) program 



at RHIC 13, £ 



It should be emphasized that energy dependence of parameters §5§ and ([5]) obtained for any centralities from most 
central (0-5%) to most peripheral (70-80%) events are similar to them shown at Fig. Q] and Fig. [5] 



IV. SUMMARY 



Two asymmetry parameters are introduced in order to investigate the initial energy evolution of possible local 
strong V '/CP violation in heavy ion collisions. Dependence on ^/snn has been obtained for parameters A a and A r 
based on the experimental correlators. The the energy dependence of absolute asymmetry for semi-central events in 
heavy ion collisions shows sharp decreasing at y / s]vjv < 19.6 GeV, almost constant behavior up to ^snn — 200 GeV 
with some decreasing at further energy increasing. Dependence iofy'sm) qualitatively corresponds to the energy 
dependence of Chern - Simons diffusion rate in Af = 4 SYM plasma at N c — > oo in strong coupling regime. The main 
features of A r energy dependence agree with behavior of A a (y/sNN) for wide energy domain under study. Thus energy 
dependencies of suggested asymmetry parameters with similar features at intermediate y/sW^j indicate possibly on 
the beginning of predominance of hadronic states versus phase of color degrees of freedom in deconfinement state in 
the matter created in heavy ion collisions with initial energies from domain ^/snn ^ 11.5 — 19.6 GeV. Therefore the 
asymmetry parameters and corresponding energy dependence seems to be useful for study of CMS and transition 
from predominance of quark-gluon phase to hadronic one. 
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FIG. 1: Energy dependence of A a for heavy ion collisions in various centrality bins: a - 20-30%, b - 30-40%, c - 40-50% and d 
- 50-60%. Experimental points are shown as following: • - Au+Au, k - Pb+Pb and ■ - U+U. Absolute asymmetry values 
are calculated based on experimental data from pi. [la. [l4T| . Smooth curves correspond to the energy dependence of Chern - 
Simons diffusion rate in the strong coupling regime for external magnetic field with strength B = (solid), B = 5T 2 (dashed) 
and B = 10T 2 (dotted). The shaded areas for T^g at B = are defined by uncertainties of T value at fixed initial energy due 
to errors of parameters in analytic function described of T(^/snn) experimental dependence fPi\ . 
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FIG. 2: Energy dependence of A r for moderate and heavy nucleus-nucleus collisions in various centrality bins: a - 20-30%, b - 
30-40%, c - 40-50% and d - 50-60%. Experimental points are shown as following: o - Cu+Cu, • - Au+Au, A - Pb+Pb and 
■ - U+U. Relative asymmetry values are calculated based on experimental data from @, 0, Q11, Q3] . 



